Non-parabolicity of the conduction band of wurtzite GaN 
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Using cyclotron resonance, we measure the effective mass, m*, of electrons in AlGaN/GaN het- 
erostructures with densities, n2D ~ 1 — 6 x lO^^cm"^. From our extensive data, we extrapolate a 
band edge mass of (0.208 ± 0.002)me. By comparing our m* data with the results of a multi-band 
k.p calculation we infer that the effect of remote bands is essential in explaining the observed con- 
duction band non-parabolicity (NP). Our calculation of polaron mass corrections - including finite 
width and screening - suggests those to be negligible. It implies that the behavior of m*(n2D) can 
be understood solely in terms of NP. Finally, using our NP and polaron corrections, we are able to 
reduce the large scatter in the published band edge mass values. 

PACS numbers: 73. 20.At, 76. 40.+b, 73. 40.-c 



The magnitude of the conduction band non- 
parabolicity (NP) in wurtzite GaN currently remains 
controversial. NP of a band can be probed by measur- 
ing the carrier effective mass m* as a function of en- 
ergy. Such experiments have been performed in the past 
in both bulk and in two-dimensional electron systems 
(2DES). The deduced band edge mass values, mg*, how- 
ever, exhibit considerable scatter. Using cyclotron reso- 
nance (CR), Drechsler et al. [ij determined mo* in bulk 
wurtzite GaN to be 0.20TOe, where rrie is the free elec- 
tron mass. Other methods such as infrared reflectivity on 
electron plasma '5] and spectroscopy on shallow donors 
d 0, in bulk GaN have yielded 0.220me < < 
0.236me. An even wider range of values for the band 
edge mass, 0.185me < TOq < 0.231me, emerge from ex- 
periments in AlGaN/GaN heterostructures. From the 
temperature dependence of Shubnikov-de Haas (SdH) os- 
cillations of 2DES, Lin et al. deduced mo* = 0.22me 
while Hang et al. reported mo* = 0.185me 0- Cyclotron 
resonance experiments on heterostructures have revealed 
0.223me < ml < 0.231me 0, This spread in mo* 
suggests that the "extrapolation" from the various ex- 
periment values m*{E) to the band edge remains poorly 
controlled. 

Using CR we have measured m* in a series of high mo- 
bihty (~20,000 cmVVsec) AlGaN/GaN structures. Our 
heterostructures are all grown by molecular beam epitaxy 
(MBE) on GaN templates prepared by hydride vapor 
phase epitaxy (HVPE). The specimens are described in 
detail elsewhere 0| . Our data cover a density range 
of 1 — 6 X 10^^ cm^^. In these 2D systems this implies 
energies from ~27meV to ~120meV above the band edge 
due to electron confinement and band filling. Therefore, 
our mass data probe the NP of the GaN conduction band 
in small steps over a wide energy range. We also perform 



extensive k.p calculations and determine that instead of 
the commonly used two-band model, a multi-band model 
is required to explain our experimental results. Addition- 
ally, our calculations on polaron correction of the effective 
mass of 2D electrons in GaN show them to be at the 1% 
level, considerably less than previously thought % A de- 
tailed comparison between our data and k.p calculation 
sets the band mass value to ml = (0.208 ±0.002)me and, 
when applied to other investigators' results, considerably 
reduces the spread in band edge mass values. 

A Fourier transform spectrometer with light pipe op- 
tics and a composite Si bolometer was used for the de- 
tection of far-infrared transmission. Magnetic field was 
applied normal to the 2D electron layer. The carrier 
density of each sample was determined in situ from the 
Shubnikov-de Haas (SdH) oscillations of the 2DES. All 
CR and SdH experiments were conducted at 4.2K. 

FigH shows the cyclotron resonance energies versus 
magnetic field, B of a sample with n2D = 2.3x10^^ cm~^ . 
All data are taken with a resolution of 0.24 meV. The in- 
set shows high field transmission spectra normalized to 
the spectrum taken at i?=OT. The solid line in Fig.l is 
a fit to the data at high fields, B>27T, and low fields, 
B<12T, resulting in an effective mass of m* — 0.228me. 
For 15T<_S<25T, there is a pronounced deviation from 
this straight line. This anomaly in the CR represents a 
recent discovery, which is being analyzed and published 
elsewhere To]. Here, we observe that this anomaly is 
limited to a finite field region outside which all CR data 
can be fit by a straight line. 

We have measured the effective mass in eleven samples 
with carrier density, n2Dj ranging from 1 — 6x10'^^ cm~^. 
In all cases, we observed either a broadening or a splitting 
of the CR line at intermediate fields but could fit the data 
away from this regime as well as the data seen in Fig. ^ 
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FIG. 1: Resonance energies vs. _B of a sample with n2D ~ 
2.3 X 10^^ cm~^. High and low field resonances can be fit 
with a single straight line, giving an effective mass of 0.228me. 
Near f 8T a level anti-crossing results in a splitting of the CR 
(see text). Inset: Transmission data for B=26, 28, and SOT, 
normalized to the spectrum at B—QT. 



FigEl shows the dependence of m* on n^ n- For compar- 
ison, we also plot data from Ref. |^ loT Tl^. The mass 
data from the different references are located in the gen- 
eral vicinity of our results but, due to their considerable 
error bar or sparsity, are difficult to extrapolate to zero 
density. The combined data of Fig|21show an increase in 
m* by ~17% as n2D changes from 1 — 9 x 10^^ cm~^. 
The rise in to* with n2D reflects the non-parabolicity of 
the conduction band of the GaN host. Simply extrapo- 
lating our closely spaced data linearly to vanishing ri2D, 
we arrive at toJ —Q.21Ame. This value is about ~8% 
lower than previously published CR data 0, • Since 
our data contain a small error bar and extend to very 
low n2D, such a simple extrapolation should already be 
quite reliable. 

Other groups have previously addressed NP in GaN. 
For example, Knap et al. explored NP using CR with 
different n2D (see our Fig|51) and accounted for the mag- 
nitude of NP using a simple two-band approximation 
9]. In such an approximation, which includes only 
coupling between the lowest conduction and highest va- 
lence bands, the effective mass varies as 



m*{E) = to;(H- 2KE/Eg 



(1) 



=1 and Eg=2,.b eV. In a 2DES the energy, E 



with K~ 

Ek + Ef, above the band minimum is composed of the 
average kinetic energy, Ek, of the electrons in the con- 
fining potential well and Ep — Trh^n2D/'n^o |l3|. the 
Fermi energy of the 2D electron system. The value of 
Ek is dependent on the form of the wave function of the 
confined electrons. In a simple triangular potential ap- 
proximation, the average kinetic energy is Ek ~ Ec/i, 
with Ec being the confinement energy of the lowest 
subband Using the more accurate Fang-Howard 
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FIG. 2; Effective mass, m*, vs. 2D electron density. Our 
data are shown as solid squares. Data from Ref. 0,1^ and 
[l^ are shown as open symbols. The dashed line is a two 
band fit to the results of Ref. |^ according to Eq^with K=l. 
The solid line represents a fit with K=2.5, which accounts for 
the infiuence of additional, higher conduction bands. Inset: 
Values of the band edge mass from Refs. Q, Q, |M |^ IB |^ 
.7. & A .12] and this work, after NP and polaron corrections. 
CR results are shown as open triangles, donor spectroscopy 
data as open hexagons and SdH data as stars. Our result of 
rriQ — (0.208 ±0.002)me is shown as a solid square. Averaging 
all data (except SdH) yields mo — 0.204 me illustrated by a 
dashed line (see text). 



variational wave function gives Ek = U^b^ / 



irriQ, where 



= ASnrnQe'^^Ndep + 5|?^2n)/e?i^ [iJl- In our samples, 
since both the MBE and the HVPE GaN are n-type, the 
depletion layer density, Ndep, can be set to zero. We use 
this two-band model with only one adjustable parameter, 
TOq , to fit the high-density data of Ref. (dashed line 
in Fig 121). Ek was computed in the triangular approxima- 
tion following the authors of Ref. |^ . Clearly, although 
this procedure can describe the original density depen- 
dence of to* of Ref. 0] due to the rather large error bars, 
it fails to account for our data. Neither a simple verti- 
cal shift of the line (a different TOq) nor the usage of the 
Fang-Howard wavefunction can resolve this discrepancy. 
What is required, is a much stronger dependence of to* 
on n2D- 

Empirically, we pursue an approach taken by Singleton 
et al. who used a modified two-band model to de- 

scribe their NP data in GaAs. The authors considered K 
in Eq^to be a second fitting parameter. The inclusion 
of a variable K > 1 into the analytic expression incor- 
porates the influence of higher conduction bands, sim- 
ulating the results of a more elaborate, multi-band k.p 
calculation Working with the Fang-Howard model 
to determine Ek we find a very good fit to our data for 
K = 2.5 and = (0.208 ±0.002)TOe. Even if we use the 
less reliable triangular approximation, the required K = 
1.9. The effectiveness of the modified expression in fitting 
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the density dependence of m* over a wide range of n2D 
demonstrates that the conduction band of GaN is more 
non-parabohc than was previously assumed 0, Q, S ■ 

Before concluding that the proposed NP model appro- 
priately describes our CR data, we need to assure our- 
selves that polaron effects are a negligible contributor to 
TO*. This mass enhancement factor results from electron- 
LO phonon coupling in polar semiconductors, such as 
GaN. A determination of a polaron mass enhancement 
in heterostructures requires inclusion of screening in 2D 
and the finite width of the electronic wave function |l7|. 
since both greatly reduce interaction between 2D carriers 
and LO phonons 0|. Following Ref. 0, we calculated 
the polaron effective mass in AlGaN/GaN heterostruc- 
tures using a a Frohlich constant, a = 0.49 0, Fang- 
Howard variational wave function, and a static Thomas- 
Fermi screening model. We find a polaron enhancement 
of less than 1% for m* for n2D = 9 x lO^^cm^^. Since 
the effect decreases with decreasing density, corrections 
for lower density specimens are smaller yet. This mass 
enhancement is considerably smaller than the 10% esti- 
mated previously for n2D = 3.1 x 10-^^cm~^ where 
screening and finite width had been neglected. A 1% 
mass enhancement due to polaronic coupling lies within 
the error bars of our CR data. 

The large spread in the value of the band edge mass 
TOq in the literature is mostly due to two reasons: the 
underestimation of the NP and the overestimation of 
polaronic corrections. Applying our NP and polaron 
corrections to the available effective mass data (Ref. 
Blli, mmilHIi 113): we reach a much more 
coherent picture for the band edge effective mass, toJ, in 
GaN, as shown in the inset to Fig. El We observe that 
the majority of the values for mg are very close to an av- 
erage TOg = 0.204TOe. However, most of the TOq data from 
SdH (displayed in stars) remain at variance from the CR, 
infra-red reflectivity and donor spectroscopy data, a fact 
that remains unexplained. 

Since m*{E) cannot be accurately represented with a 
two-band model, a five-band k.p calculation in the 
zincblende approximation was performed to model our 
data. Figure 13 shows the results. The material param- 
eters employed using the Koster notation are the spin- 
orbit splitting, Ag, of the Fs conduction band and the 
momentum matrix elements, and X^P^, coupling the 
Fi conduction band with the valence and conduc- 
tion bands. They are not precisely known. The formal- 
ism proposed by Carlos et al. jjjj and used by Bayerl 
et al. |2i| in a 5-band model to relate the parameters 
was utilized to calculate m*{E). However, as seen in 
FigEl the parameters chosen by Bayerl et al. lead to a 
non-parabolicity even lower than in a two-level model. 
We find that the experimental data points can only be 
matched if an extra parameter C, taking remote bands 
into account at k=Q, is included [i^. If one assumes 
C= -1.5 (compare with C=-2 for GaAs 0), is in 
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FIG. 3: Variation of m* with energy according to k.p calcla- 
tions, assuming m3=0.208me. The two-band model and the 
five-band k.p results of Bayerl et al., which neglect remote 
bands both underestimate the energy-dependence of m*. The 
experimental data (solid squares) can be explained by multi- 
band k.p calculations that include the influence of remote 
bands (dashed lines). 

the range 25-27 eV with «0.33-0.48 and A[, «120-180 
meV. This set of parameters is close to the "best" set of 
values chosen by Kennedy et al. |2^ and confirm that 
p2 is relatively large for GaN (P^ ^26 eV). 

In conclusion, our CR experiments set the conduction 
band edge mass in GaN to TOg = (0.208±0.002)TOe. Using 
our determination of the nonparabolicity and reviewing 
the polaron correction we reach much better agreement 
between several published data for TOq . 
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